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Abstract 

We consider the supersymmetric approach to gaussian disordered systems like 
the random bond Ising model and Dirac model with random mass and random 
potential. These models appeared in particular in the study of the integer quantum 
Hall transition. The supersymmetric approach reveals an osp(2/2)i affine symmetry 
at the pure critical point. A similar symmetry should hold at other fixed points. 
We apply methods of conformal field theory to determine the conformal weights 
at all levels. These weights can generically be negative because of non-unitarity. 
Constraints such as locality allow us to quantize the level k and the conformal 
dimensions. This provides a class of (possibly disordered) critical points in two 
spatial dimensions. Solving the Knizhnik-Zamolodchikov equations we obtain a set 
of four-point functions which exhibit a logarithmic dependence. These functions 
are related to logarithmic operators. We show how all such features have a natural 
setting in the superalgebra approach as long as gaussian disorder is concerned. 
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1 Introduction 



There is growing evidence that some disordered systems at criticality share two unusual 
features: The existence of logarithmic operators in the spectrum of the theory and the 
existence of an infinite number of relevant operators with negative conformal dimensions. 
Logarithmic operators seem to be connected to hidden continuous symmetries. In Jl|] 
logarithmic operators were found to generate a change in a coupling constant of the 
effective WZNW action for SU(2r), r — > 0, obtained after use of the replica trick. The 
existence of an infinite number of conformal operators |2| means that there is an infinite 
number of relevant perturbations which render the critical point unstable. Such a set of 
conformal dimensions is also related to the phenomenon of multifractality. 

Our aim in this paper is to show that both features follow naturally from a supersym- 
metric treatment of the disorder. We consider the two-dimensional random bond Ising 
model at criticality and the random Dirac model in 2 + 1 dimensions. These gaussian 
models allow a supersymmetric formulation of averages over disorder ||. One identifies 
a global superalgebra symmetry of the effective action and assumes it is enhanced to an 
affine symmetry at a new critical point || |J, as happened at the pure critical point. One 
then derives the Sugawara stress-energy tensor and obtains the set of conformal dimen- 
sions associated with the primary fields. Additional constraints such as locality allow to 
further restrict the operator content of the theory. This provides a class of, possibly disor- 
dered, critical points in two spatial dimensions. The specific structure of the superalgebra 
implies that such dimensions can generically be negative. The presence of logarithmic 
operators is also straightforward from the point of view of this algebra. Unlike ordinary 
Lie algebras, superalgebras have indecomposable (not fully reducible) representations. It 
is then possible to show on general grounds that such representations imply the existence 
of logarithmic operators and logarithms in correlation functions. 

This paper is organized as follows. In section 2 we briefly recall the supersymmetric 
approach for the random bond Ising model and the random Dirac model. In section 3 we 
write down the osp(2/2)k current algebra and find the Sugawara stress-tensor. In section 
4 we briefly review the representations of osp(2/2) and obtain the conformal weights of 
primary fields. In section 5 we derive the quantization of k from locality constraints. In 
section 6 we show how indecomposable representations lead to logarithmic operators and 
consider an atypical representation. In section 7 we obtain some logarithm-dependent 
four-point functions. We conclude in section 8. 

2 The supersymmetric method for disordered sys- 
tems 

The supersymmetric method is applicable to models which are gaussian at fixed disorder, 
but it provides a good starting point for disentangling properties of a large class of disor- 
dered conformal field theories. The study of the random bond Ising model reveals that an 
appropriate algebraic framework for studying gaussian disordered systems at criticality 
should be based on affine Lie superalgebra with zero superdimension. In the case of the 
random Ising model this algebra is osp(2N/2N). The fact that the algebra has an equal 
number of bosonic and fermionic generators ensures that the Virasoro central charge of 
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the Sugawara stress-tensor vanishes. This is needed by construction for a disordered sys- 
tem. The vanishing of the central charge does not imply the triviality of the theory since 
it is not unitary. The relevance of affine Lie superalgebras was independently realized in 
refs. [|, 1- 

We now briefly recall the supersymmetric method for the study of the random bond 
Ising model and the random Dirac model. These models have been analyzed using the 
replica method in refs. [|], ||, 0. The random bond Ising model has also been studied in 
the context of massless scattering theories 0. 

The random bond Ising model 

In the scaling limit, near criticality, the Ising model is described by a massive real 
Majorana fermion with mass m ~ r = ^T 7 ^ were Tq is the critical temperature. In the 
presence of disorder, i.e. when the coupling constants between sites belong to a random 
set, the mass becomes a function of space position. The random Ising model is defined 
by the action (z — x + iy): 

cpx 

im(x)ibib) (V, 



/(X X / 
— - [ipdzip + ipd z ip + 
4:71 \ 



where ip and ip are grassmanian fields. The mass m(x) is chosen to be a quenched random 
variable with a gaussian measure: 

P[m] = exp (-^- J ^-(m(x) - mf^j . (2) 

The energy operator e(x) = itp(x)ip(x) has dimension one. The Harris criterion tells 
us that randomness in the bond interaction is marginal in the 2d Ising model. It turns out 
that it is not exactly marginal but only marginally irrelevant. At criticality the disorder 
only induces logarithmic corrections to the pure system. 

In order to compute averages of products of correlation functions one introduces a 
number of copies of fermions and of their supersymmetric partners equal to the number 
of correlation functions in the product j3], |J. One then rewrites these averages as fermionic 
and bosonic path integrals. For the Ising model with m = one obtains the following 
effective action for the disorder average of the product of two correlation functions: 

d^x d^x 
= S. + f / — $ pcrt (3) 

with 



$pert = -^M>+ +V1-V1) ■ (4) 



The ip are complex fermions and rj and 7 are complex bosonic fields. This action can 
be viewed as a perturbation of the (non-unitary) conformal field theory specified by the 
action S*. This fixes the normalization of the fields to be: 

if>-(z)if)+(w) — , j(z)r)(w) ~ — - — . (5) 

z — w z — w 
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The central charge of the Virasoro algebra is zero. Note that since the fermions ip± have 
dimension one half the perturbing field $ per t has dimension two. It is therefore marginal. 

The conformal field theory specified by is invariant under an affine supersymmetric 
algebra whose conserved currents are: 

G±(z) = r](z)ip±(z) , G±(z) = 7(z)V>±(z) , 
K(z)=r, 2 (z) , K(z)= 1 2 (z), 
J(z)=:^_(^+(z): , H(z)=:-y{z) V {z): . (6) 

The dots refer to fermionic and bosonic normal ordering. There are four fermionic cur- 
rents, G± and G±, which are generators of supersymmetric transformations, and four 



bosonic ones. They form a representation of the affine osp(2/2) |12], [13] current algebra 
at level one. 

The perturbing field can also be written as a bilinear in the currents, 



$pcrt 



JJ-HH + -(KK + KK) + G_G+ - G_G+ + G+G_ - G+G_ 



(7) 



In other words the perturbation (|j) is a current-current perturbation. Therefore the 
action S eff preserves a global osp(2/2) symmetry. 

The random Dirac model 

The random model of Dirac fermions has been introduced in connection with the 
quantum Hall transition ||. Its action is: 

d 2 



— ty_d s lf> + + l/>_d g l{>- 



- 1>-1>+) + *^(^+ + MS) ■ (8) 

The random variables m and V have a gaussian distribution with widths Qm and gy. 
We denote by $m and $y the perturbing fields coupled to the constants gu and gy after 
averaging over the disorder. In the two-copy sector we can write them in terms of the 
osp(2/2) currents: 

$ v = 2HH — 2 J J + TCK + KK + 2G_G + — 2G + G_ + 2G_G + - 2G + G_ , (9) 
$ M = 2JJ-2HH + KK + KK + 2G_G + -2G + G_-26_G + +2G + G^ . (10) 

We see that $ M = ^pert- 

Both <3> M and $y preserve a global osp(2/2) symmetry, whose generators are 



H + H , Jo+_Jo, K + K , K +J{ 

G + o + G + o G + q + G + o , G_o + G_o , G_o + ^-o 



for the perturbation $Mi and 



H + H , J +_J , K -K , K -_K 

G+o + G + o j G + o — G + o j G_o — G_o > G_o + G_o 



(12) 
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for the perturbation <3>y. Here the index denotes the zero modes of the currents 
J a (z) = Y!m Jn z ~ n ~' 1 ■ ^ both perturbations are present, the interaction preserves an 
u(l/l) symmetry generated by 

Ho + Ho, Jq + Jq, G + q + G + o , G-o + G_o • (13) 

The perturbative study of the Dirac theory with a random potential and a random 
mass is very similar to the perturbative study of the random mass Ising model. However 
the crucial difference is that contrary to the randomness of the mass, the randomness of 
the potential is marginally relevant. At gm = 0, the one-loop beta function is given by 
9v = $>9v- This means that gy grows at large distances. The infrared fixed point has not 
yet been determined 0. Should the affine symmetry, present at the UV fixed point of 
the pure system, be restored at the IR point, then we believe this point should belong to 
the set we find in section |[ 



3 The current algebra approach 

The critical free theory described by the action has a current algebra symmetry, 
osp(2/2)i, on both its holomorphic and antiholomorphic sectors. When randomness is 
introduced in the pure system the effective action is no longer critical, and the current 
algebra symmetry is reduced to a global osp(2/2) symmetry. However at a new critical 
point the conformal invariance restores the current algebra symmetry. The value of the 
level k is not preserved by the renormalization flow. Therefore its value at the IR fixed 
point could be different from its value at the UV point. It is possible to extract properties 
at the new critical point by studying the current algebra and its associated stress-energy 
tensor. In this context we consider osp(2/2)k at arbitrary k and obtain the Sugawara 
tensor. 

3.1 The osp(2/2) k algebra 

In this section we write down the singular terms of the Operator Product Expansions [|Kj] 
satisfied by the currents of the affine osp(2/2) algebra at level k. The non-trivial OPE 
of the currents are easily obtained from the currents (|6|). The level k appearing below is 
equal to 1 for these currents. We find: 

J( Z )J( W )~-J^— ; H(z)H(w) 



(z — w) 2 > (z — w) 2 

J(z)G±(w) ~ -^—G±(w) ; J(z)G ± (w) ~ -^—G ± (w) 
z — w z — w 

H(z)G±(w) —G±(w) ; H(z)G±(w) ~ -^—G±(w) (14) 

z — w z — w 

H(z)K(w) — K(w) ; H(z)K(w) — 

z — w z — w 

G ± (z)G T (w) ~ —*+-!—(H(w)±J(w)) 
{z — w) z z — w 
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G-(z)G+(w) 
K(z)G±(w) - 



K{z)K(w) 
-K{w) 



2k 



(z — w)' 



+ 



z — w 



-H(w) 



z — w 
-2 



G±(w) 



G4z)G + (w) 
K(z)G±(w) 



z — w 
2 



-K{w) 



G±(w) 



z — w z — w 

Requiring the algebra to be associative constrains the possible central terms extensions 
to the ones appearing in the OPE's (14). We can rewrite these OPE's in a more compact 
form as 



,.ab 



J a {z)J b (w) 



z - w 



+ f 



ab 



J c {w) 



z — w 



(15) 



The f c are the structure constants of the Lie superalgebra osp(2/2), and k is propor- 
tional to its non-degenerate Killing form. 



3.2 The Sugawara stress-energy tensor 

We now construct the Sugawara stress-energy tensor |TT|. It is bilinear in the currents 
defined in the previous section. Because of singularities which appear at coinciding points 
one has to consider a regularized version where normal ordered products of currents ap- 
pear. This is equivalent to a point splitting procedure where the singular parts appearing 
in the OPE's of the currents are subtracted. We take the usual definition for the normal 
ordered product of two fields A(z) and B{w): 

:AB:(w)=<f (16) 
Jw 2iri z — w 

It should be noted that this ordering prescription does not coincide in general with the 
Wick prescription used in (|). The Sugawara stress-energy tensor T(z) is given by 

T(z) = -K ab : J a {z)j\z) : (17) 

where is the inverse of K ab . The constant k is determined by requiring J a (z) to be a 
primary field of conformal weight one: 

T(m „ J^) + f!£!M . (18) 

[z — w) 2 z — w 

When calculating OPE's one has to add a minus sign each time two odd generators are 
permuted. We find 

T(z) = J^((-H(z)H(z)-J(z)J(z)-~(K(z)K(z) + K(z)K(z)) 

+ G + (z)G^(z) -G-(z)G+(z) + G-(z)G+(z) -G+(z)G-{z) :) . (19) 
The Sugawara tensor also satisfies 

T(z)T(w) ~ 1 + 2 1 ; + . 20 

(z — (z — wy z — w 
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For a general Lie superalgebra, it is easy to show that the Virasoro central charge c is 
proportional to the super dimension; this is the difference between the number of even and 
odd generators. Therefore c vanishes for the osp(2/2) algebra, and more generally for any 
superalgebra with an equal number of even and odd generators. More precisely one has: 

c = 2k^ = 2k S -^ = 0. (21) 

K K 

Because c vanishes T(z) is a primary field, contrary to the case c ^ where T is just a 
level-two descendant of the unit operator. 

4 Primary Fields 

We briefly describe the representations of osp(2/2) and find the corresponding conformal 
weights. 

4.1 Some osp(2/2) representations 

Unlike ordinary Lie algebras, the are two types of representations for most superalgebras. 
The typical representations are irreducible and are similar to those of ordinary Lie algebras. 
The atypical representations have no counterpart in the ordinary Lie algebra setting. They 
can be irreducible or not fully reducible (read reducible but indecomposable). 

The superalgebra osp(2/2) is isomorphic to the superalgebra spl (2/1). The repre- 
sentation theory of the latter algebra was studied in fl2|, [13| . The quadratic Casimir of 
osp(2/2) is 

C 2 = - (7f 2 - -(KK + KK) - J 2 + - + - G+G^j . (22) 

The four even generators K, K, H, J form a su(2) ®u(l) subalgebra. The correspondence 
with the notation of p] is: Q + = -K/2, Q_ = K/2, Q 3 = -H/2, B = J/2, V+ = 
G+/V2, V_ = -G+/V2, W + = -G_/V2, W_ = G^/^[2. 

Let b and ^3 be the eigenvalues of B and Q3. Generically, a representation (b,q), 
b E C , q = 0, |, 1, |, contains four su(2) © u(l) multiplets: 

\b,q,qs) , 13 = ~q, -q + l,...,q- l,q if q > , (23) 
.,11. 1 311 

' h+ 2' q ~2' qz ' ' q3 = ~ q+ 2 , "' ,q ~ 2 ,q ~ 2 g ~2' ^ ' 

,,11. 1 311 

' b ~2 ,q ~2 ,q3 ' ' q3 = ~ q+ 2'"' ,q ~ 2 ,q ~ 2 q ~2' ^ 
\b,q-l,q 3 ) , q 3 = -q + l,...,q-2,q-l if q>l. (26) 

The action of the four even generators on these multiplets is the one implied by the 
notation. The four odd generators mix the different multiplets. The vector \v) = \b,q,q) 
is a highest weight vector, i.e. it satisfies 

K\v) = G+\v) = G-\v) = . (27) 
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The quadratic Casimir of this representation is C2 = 2(g 2 — b 2 ). 

If b ^ ±q the representation is denoted by [b, q] and is typical; the quadratic and 
cubic Casimirs do not vanish. All the vectors in the representation can be obtained 
from the highest weight vector \v) by applying on it polynomials in the generators. The 
representation [b,q] has dimension 8q. The representation [0, 1/2] is four-dimensional and 
contains one spin 1/2 and two spin multiplets. 

Atypical representations 

When b = ±q several kinds of atypical representations arise. Both Casimirs vanish, 
and yet these representations are not the trivial one-dimensional representation. One kind 
has dimension 4q + 1 and is denoted by [q]±. To obtain [q] + (resp. [q]-) one drops the 
two multiplets (|25|) and (|26|) (resp. (|24]) and (|26|)). These representations are irreducible. 

Atypical indecomposable representations 

Generally, they are semidirect sums of atypical irreducible representations. They can 
contain two, three or four terms and they arise in tensor products of irreducible represen- 
tations. An interesting example, containing four irreducible representations, arises in the 
tensor product of two representations [0, 1/2]. The result is the direct sum of [0, 1] and 
an eight- dimensional representation which is the semidirect sum of [1/2] _, [1/2] + and two 
[0] representations. We call it [0,-1/2,1/2,0] (see fig. 1 b). The vector s is invariant; 
it is annihilated by all the generators. The quadratic Casimir Ci vanishes on all states 
except t: C 2 t = 4s. We have 

5 = o> - |o, |> |0, -|> + [0, -i> ^) [0, |> + 1 - |, 0)®|i,0) 

* = I 0> <S) I — |, 0> -h |0, |> <2) |0, — i> — [0, — i> ^ | , i> -h I — 0> <2) |i, 0>) 

4.2 Tensor product of osp(2/2) representations 

The tensor product of two irreducible representations of a superalgebra is not necessarily 
completely reducible. Ref. [|ll| gives a sufficient condition for a tensor product of two 
osp{2/2) representations to be completely reducible, and the irreducible components. The 
result is the following 

[b, q] <g> [b', q'] = [b + b', q + q 1 ] © [b + b', q + q' - 1] © . . . © [b + b' , \q - q'\) 
©[& + b', q + c( - 1] . . . © [b + b', \q - q'\ + 1] 
®[b + b' + 1/2, q + q'- 1/2] . . . © [b + b' + 1/2, \q - q'\ + 1/2] 
®[b + b' - 1/2, q + q'- 1/2] . . . © [b + b' - 1/2, \q - q'\ + 1/2] (28) 
if ± b > q > f , ±b' > q' > \ . 

For the other values of b, b', this expression still gives the correct content of su(2) © u(l) 
charges of the tensor product. This expression also permits to obtain some information on 
the reducibility: if two components have different Casimirs, then they belong to different 
(maybe indecomposable) representations. Vanishing Casimirs are, generally, signs of some 
pathologies. 
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A b 





a) 



b) 



Fig.l: a) The tensor product of the representations [b, 1/2] and [b', 1/2]. If b > 1/2, b' > 
1/2, the result is the direct sum of [b + b' , 1] (white dots), [b + b' + 1/2, 1/2] (black dots) and 
[b + b' — 1/2,1/2] (grey dots). The arrows represent the action of the generators, b) The 
representation [0, —1/2, 1/2, 0]. The vector t (black dot) is cyclic and s (white dot) is invariant. 
The action of the fermionic generators G+, G_, G+ and G_ is represented by arrows. 



4.3 osp(2N/2N) representations 

When calculating the disorder average of product of correlation functions for the Ising 
and Dirac models, we can introduce an arbitrary number of copies of fermions and bosons. 
One is led is to consider osp(2N/2N) as symmetry. The theory is consistent if the results 
of the calculations do not depend on N. One step in this direction is to show that the 
conformal dimensions of the primary fields are independent of N. 

The main difference between the superalgebras osp(2/2) and osp(2N/2N) with N > 1 
is that the even subalgebra of the first one is not simple. This is why representations 
of the former are indexed by a continuous parameter, while the representation of the 
latter are characterized only by discrete variables. According to the classification of the 
superalgebras |14[], osp(2/2) is the first of the series C(N + 1), while osp(2N/2N) is a 
superalgebra of the type D(N, N). 

The finite dimensional typical representation of the superalgebras were classified by 



Kac [[14 1 . For D(N,N), the typical representations are characterized by their highest 
weights, A = Y^iLi a i$i + Hf=i a N+i^i ; <^ and q , i = 1, N form an orthogonal basis with 

0. The 'numerical marks' a, satisfy the following 



conditions: 

i) Oi G Z+, i N, 

ii) 3 = a N- ajv+i - ••• 

iii) a N+ j +1 = ... = a 2 N 



a 2N-2 — 2( fl 2iV-l + a 27v) £ Z +j 

0, if j < N - 2; a 2N -i = a-2N, if j 



N-l. 
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The first two conditions express the fact that A is a dominant weight for the even 
algebra Sp(2N) <g> 0{2N). 

The value of the quadratic Casimir for the representation with highest weight A is 
given, up to a normalization, by (A, A + 2p), with p = Po — p\ and po(i) is half the 
sum of the even (odd) positive roots. Considering the positive roots associated to the 
distinguished Dynkin diagram (the Dynkin diagram with a single fermionic root), we 
obtain 

N 

P = ^2( N - 0( e i - &V--i+l) • 
t=l 

We normalize the value of the quadratic Casimir of the adjoint representation (2, 0, 0) 
to be 1; the Casimir of the representation (ai, a 2 N) will be 

1 N 

C 2 {K) = -Y,[{a l -i + lf-{a N+i + N-i) 2 . 

4 i=l 

For the fundamental representation (1, 0, 0) we have C 2 — 1/4, independent of N. 

In order to obtain compatibility between the results obtained in the osp(2/2) and 
osp(2N/2N) frameworks, the continous parameter b of the osp(2/2) representations has to 
be constrained to the discrete values compatible with the value of the discrete parameters 
of the osp(2N/2N) representations. 



4.4 Conformal dimensions of the primary fields 

The primary fields <j) a (w) are highest weight vectors of affine osp(2/2)k and of the Virasoro 
algebra. In terms of operator products one has: 



z — w 



(29) 



T{z)(p (w) ~ A- - + r . 30 

[z — w) 2 (z — w) 

Using ( JT5| ) and one then shows that the matrices T a form a representation of the 
osp(2/2) algebra. But this applies to any algebra. Using the Sugawara form (fl7f) of 



the stress-energy tensor and (|T^) one obtains the conformal weights corresponding to a 
particular representation. For the representations [b, q] and [q]± the Casimir C 2 is diagonal 
and we obtain: 

A(M = U b T«T» = 2 -^l . (31) 

We are also able to calculate the conformal dimensions of the primary fields for all N. 
Due to the fact that the dual Coxeter number of osp(2N/2N) equals 1 for all N, these 
conformal dimensions do not change if we change the number of copies to N' > N; for 
example 

. _ 1 

A (1 , ,...,o) - j^k 

is the conformal dimension of the fields in the fundamental representation. 

The conformal dimensions are potentially negative; this reflects the non-unitarity of 
the theory we considered. 
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5 Local Operators and Quantization of k 



As seen in section two the interactions generated by the random mass and random scalar 
potential conserve diagonal osp(2/2), J£ = Jq ± Jq. We expect to find in the spectrum 
of the theory, at the new fixed point for the potential interaction, particles which carry 
the same charges H x and J x as the original ones (bosons and fermions). The fields 
corresponding to these particles are no more holomorphic or antiholomorphic, but they 
can be made up of the new primary fields. Suppose that they are of the form 

<j)(z, z) = [6)?] (z)(j) [Vjq/] (z) (32) 

where 4>[b^(z) (4>[b' ,q']{z)) are primary fields belonging to the representation [b,q], respec- 
tively [b',q'] of the two chiral algebras J a (z) and J a (z). 

The charges of the bosons and fermions with respect to the currents J x are: 





Jr 


H x 




1/2 





iM*) 


-1/2 





tO) 





1/2 


7]{z) 





-1/2 



We have to identify the representation [0, 1/2] in the tensor product of the representations 
[b,q] and [b',q']. Using equation (p8| ) we obtain the following possibilities 

a) q = q>, b + b' ±1/2 = 0, b) \q - q'\ = 1/2, b + b' = . (34) 

As in 0, in order to obtain a sensible theory, we impose the constraint of locality on 
the fields (f>(z, z). The OPE of two such composite fields can be written as 

O a (z, z)O b (0, 0) = X; C c ab z A «- A «- A * z A «- A "- A »O c (0, 0) . (35) 

c 

The operators O a and Ob are said to be mutually local if the OPE is single-valued, which 
means that the spin difference is an integer 

s c - s a - s b = A c - A a - A 6 - A c + A a + A 6 e Z . (36) 

Let us take O a (z, z) = Of,(z, z) = ^^{z)^ ^{z) . We evaluate the left and right chirality 
OPE's separately. The operators appearing in the OPE for the left chirality are the 
operators associated to the representations a appearing in the tensor product [6, q) ® [b, q] 
(see eq. @) 

~ E^ A ^ 2A[M ^(°) 

a 

and similarly for the right chirality, with [b, q] replaced by [b', q'\. We ignored the possible 
existence of logarithms in this OPE (when two conformal dimensions are degenerate). 
We also supposed that all the primary fields allowed by the osp(2/2) selection rules (and 
their descendants) appear in the OPE. For a Kac-Moody algebra this is not always the 
case; the operator content of the theory depends on the value of the level k. In fact the 
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calculation shows that the constraint of locality depends very weakly on how many of the 
operators (fi a really appear in the OPE. 

The condition s c — s a — s b G Z is fulfilled for all the possible operators O c (0, 0) in fl3"5j ) 

if 

2l(4q - 1) G Z , 2/(4g - 2) G Z , Al{q - b) G Z , (37) 
2l(4q' - 1) G Z , 2Z(4g' - 2) G Z , Al(q' - V) G Z , (38) 
and 4l{q 2 - b 2 - q' 2 + b' 2 ) G Z , (39) 

where I = 1/(2 — k). With the two possibilities in (^) we obtain the following restriction 
on k: 

k = 2 - - , I G Z . (40) 



6 Logarithmic operators 

We argue that the appearance of logarithmic operators is a general feature of theories 
with an underlying superalgebra symmetry, as long as indecomposable representations 
are concerned. We give as an example the occurrence of logarithmic correlations for 
osp(2/2). 

The first logarithmic correlation functions in conformal field theory were obtained by 
Rozansky and Saleur in Jl5| in the framework of the GL(1, 1) WZW model. Later Gurarie 
[L6| pointed out that the logarithms are related to non-diagonalizable Virasoro generators. 



See also |1J for disordered models, [|17]], and |T8| in relation to the fractional quantum Hall 
effect. Suppose we are given two operators Ai(z), A 2 (z) with: 

T(z)AM ~ ^4 l( 0) + ^, (41) 

T(z)A 2 (0) ~ ^A 2 (0) + ^^ + 4^i(0). (42) 

z z z Z A 

Then they transform under the dilatations z' = a z as 

A'^z') = a-^A^z) , A' 2 (z') = a- A A 2 {z) + alogaA 1 {z) . (43) 

And more generally one has: 

[L n , A 1 (z)} = {z n+1 d z + A(n + 1)^)^) , (44) 
[L n , A 2 (z)} = {z n+1 d z + A(n + l)z n )A 2 {z) + a{n + l)z n A x {z) . (45) 

This gives rise to differential equations for the correlation functions of the operators A\ (z) 
and A 2 (z), through the conformal Ward identities for n = —1,0,1. Let us determine 
the two-point functions of these operators. The n = differential equations for these 
correlation functions are 



(zd z + 2A)(A 1 (z)A 1 {0)) = , (zd z + 2A)(A 2 (z)A 1 {0)) + a(A 1 (z)A 1 (0)) = , 
(zd z + 2A)(A 2 (z)A 2 (0)) + a(A 1 (z)A 2 (0)) + a(A 2 (z)A 1 (0)) = , 
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while for n = 1 we have 



This implies 



(zd z + 2A)(A 2 (z)A 1 (0)) + 2a(A 1 (z)A 1 (0)) 
(zd z + 2A)(A 2 (z)A 2 (0)) + 2a{A 1 (z)A 2 (0)) 



(A 1 (z)A 1 (0)) = 0, (A 2 (z)A 1 {0)) = (A 1 (z)A 2 (0)) 
(A 2 (z)A 2 (0)) = {-2aC log z + C')z~ 2A 



0. 

. 



Cz 



-2A 



with C, C arbitrary constants. This system of correlation functions appeared in [[U| 
||. Thus if the operators A u A 2 satisfy the OPE (|T]-||) they can have logarithmic 
correlators, and A 2 (z) is the logarithmic operator related to Ai (z). 

In the following we give an example of such operators appearing in the indecompos- 
able representations of the osp(2/2) superalgebra. Consider fields 4>i(z) in the atypical 
indecomposable representation [0, —1/2, 1/2,0] (fig. lb). The structure of the generators 
T a is 

A A' s t 
X ( t a * ^ 
t' a * 
* * 

V o o ooy 



rpa 



X' 

s 

t 



(46) 



where A = [1/2] + , A' = [1/2] _ and the stars stand for possible non-zero matrix elements. 
The action of the currents on the fields takes the form 



(o) 



r(z)Mo) = T^ + T^ 



(47) 



where the index 1 refers to the invariant subspace. This form is typical of indecomposable 
representations of any algebra. We find the most singular term arising from the J — <fi 
contractions in ( J a (z) J b (w)4> 2 (0) Y[ O). Using the definition 



T(z) 



lim — 

z-*w K 



K nh J a (z)J b (w) 



K a i,K 



ah 



(z — wy 



we conclude that 

T(z)M0) 
T(z)M0) 



4<Mo) + -&M0), 

z z z 

A i K^n + T-Tt) 
-^<MUJ H 

Z A K Z z 



>i(0) + -90a(O) 

z 



(48) 

(49) 
(50) 



where A = ±K ab T?T?. This non-diagonal action of the Virasoro algebra is the direct 
consequence of the existence of an indecomposable representation at the affine level. For 
the representation at hand one gets C 2 1 ~ s, C 2 s = 0, and for the fields <p s and <p t with 
zero osp{2/2) charges: 



T(z)<j) 8 {0) , T(z)<p t (p) ~ a — — H 



(51) 
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According to the previous discussion their correlation functions are 



(^(*)&(0)> = , (Mz)M0)) = C , (&(z)&(0)) = ~2aC log z + a . 

We can also solve the differential equations for (<j> s ) and ((f>t) to find (4> s ) = 0, (<f> t ) = c. 

The indecomposable representation we considered appear in the tensor product of two 
fundamental representations. Thus logarithmic terms are expected to appear in the four 
point functions (if the logarithmic operators are not identically zero, as seems to be the 
case for the free bosons/fermions k — 1). 



7 Four-point functions 



We now calculate a set of four-point correlation functions for primary fields in the [0 
representation of osp(2/2). As explained in section 
dependence. 



' 2i 



we expect to find a logarithmic 



7.1 Differential equations 

Consider the basis 

|i> = |o,f>, |2>H- 



1,0), |3) = |i 0), |4) = |0,-|) 



(52) 



where |1) and |4) are even (bosonic) and |2) and |3) are odd (fermionic). These vectors 
are labeled by the eigenvalues of B and Q3. The odd generators are given by: 



V + = V2 



V2 



( 


e 





^ 






( 

























, w. 


= V2 


-p 




















a 




















V 








J 






\ 








7 


J 


f 





7 


^ 




1 


f 











°\ 













, V- 


= V2 
































-a 











V 








J 




\ 


\ 





€ 





/ 



(53) 



(54) 



The four parameters appearing in (53-54) are constrained by 

4«7 = 1 , 4f3e = 1 . 



(55) 



Thus there are two free parameters which correspond to arbitrary relative normalizations 
of the su(2) doublet (|1), |4)), and the two singlets \2) and |3). 

The vacuum invariance under the zero modes of the affine algebra, or equivalently 
invariance under global gauge transformations, imply the following linear equations 



(56) 



i=l 



It is understood that minus signs arise each time we permute two fermionic objects, and 
omit the antiholomorphic dependence and indices of the fields. Similar equations hold for 
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the 'right' (antiholomorphic) chiral generators. These equations are solved in terms of 
invariant tensors for the representations at hand. We have found the explicit expressions 
for three linearly independent tensors which we give in the appendix. 

Requiring the fields of a specific module to be both affine and Virasoro primary gives 
a set of equations relating the fields. Such equations give partial differential equations, 



the Knizhnik-Zamolodchikov equations p9[, for the correlation functions. One gets 



i_d_ 

2dz^ 



E 



rparpb 



K j=2 Zl z i 



>i(2i)n<M%)>=°- 

i=2 



(57) 



After a conformal transformation this equation reduces to a matrix differential equation 
for the correlators 



z) = lim w 2A4 (0 il (^)0 t2 (O)0 l3 (l)0 i4 («;)) . 



(58) 



One gets 



where 



1 



2 

^12^34 



dF(z) 



i /i 



K \Z 



-V 



Q)F(z) 



(59) 



(60) 



-rv r ~ r \ _ -2Ai -A1-A2-A3+A.A1-A2+A3-A.A1+A2-A3-A4 cv r N 
^-^1, ^2, Z3, Z4) — Z u Z 23 Z 24 Z u r \Zj . 



^23^41 

Similar equations hold for the antiholomorphic sector. The complete correlator is given 
by 

y -2Ai„ -A1-A2-A3+A4 A1-A2+A3-A, 
'14 z 23 ^24 z 34 

For the foregoing representations Aj = 1/ (2x) and the Zij prefactor reduces to z\± z 2 £ 
where 

x = 2-k. 

The correlators have the osp(2/2) x osp(2/2) invariant decomposition 



^plH2«3«4JlJ2J3J4/^ g\ 



*A 1 B t<AB{Z,Z) 



(62) 



A,B=1 



The tensors are given in the appendix and the nine scalar functions satisfy the differential 
equations: 



x- 



OF 

dz 



z 



1 



1 



Q 



(63) 



where F now denotes the vector (F±b, F2B, F^b) for all B. There are similar equations 
for the antiholomorphic dependence, for the vector (Fai, Fa2, Fas)- We give the matrices 
P and Q in the appendix. Let /a(z) = Fab(z) for any given B. It is straightforward to 
reduce this first-order matrix differential equation to the following set of equations:[] 

x 3 z 3 {l - zfd 3 f 3 {z) + x 2 {l + 2x)z 2 {l - zf(l - 2z)d 2 h(z) + xz{\ - z)[-l - x 
+2xz - 2x(2 + x)z{l - z)]df 3 (z) + {-x - 1 + 2z + 2xz{\ - z))f 3 {z) = (64) 



1 It is amusing to note that the third-order differential equations for fi(z) and f%(z) reveal two apparent 
singularities at z = 2 and z = i. 
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Hz) = - 4eixz{1 _ z) {* 2 D 2 h{z) + 2x(l - z)Df 3 (z) + (1 - 2z)f 3 (z)) (65) 

Hz) = -^-(xDf 3 (z)-f 3 (z)) + (z-2)f 2 (z) (66) 
4e7 

where .D = z(l — z)d. As expected the equation for f 3 is fuchsian with three regu- 
lar singularities at 0, 1 and oo.Q The corresponding indices are (— x -1 ,— x _1 ,l + x -1 ), 
(— 1 — x _1 ,x _1 ) and (0,0,1 + 2x~ 1 ). The degeneracy of some indices, i.e. when 
the difference of two indices belongs to Z, signals the eventual presence of logarithmic 
solutions. Thus one is naturally lead to distinguish between different values of x's. 

We were able to express the functions /$ in terms of generalized hypergeometric func- 
tions. We obtained: 

F^{z) = h{z) = \{l-xD)F-{z), 
x 3 z 3 (l - zf&F- + x 2 (3x + l)z 2 (l - z) 2 {\ - 2z)d 2 F~ 
+z(l - z)(x 2 (x + 1)(1 - 6z(l - z)) - x)dF~ - (1 - 2z)F~ = . (67) 



Equation ( p7|) is fuchsian with three regular singular points at 0, 1 and oo. The indices 
at and 1 are (— — -, -). The solutions of equation (|5"7| ) can be written as generalized 
hypergeometric functions: 

F -(z) = (zil-z))' 1 * 3 F 2 (l, --,1 --;1,--; Az(l - z)) (68) 

\2 x x x J 

Fr(z) = F -(z)log[iz(l-z)} + (z(l-z))- l *J2 c n(*z(l-z)) n (69) 

n>0 

/I 1 1 2 2 2 A 
F 2 -^ = 2 1-2))= 3^2 -1 + -- + -;1 + -,1 + -; Az(l-z)) 70 

\x x 2 x x x J 



where 



T? ( U U \ — ( a l)n( a 2)n(03)n n / \ / , i \ / , 

3 F 2 (ai,a 2 ,a 3 ;6i,6 2 ;z) = > TTTTT^ — i - 2 > {a) n = a(a + 1) ■ ■ ■ (a + n 

(6i) n (& 2 )„n! 

T(n + s + |)I> + g - |)T(n + 5 + 1- |)T(s + lffls + l~j)\ 



r(n + s + 1) 2 I> + s + 1 - 2)r( a + i)r(s - ±)r(s 1 1 1 



s=0 



It must be noted that although the solutions look symmetric under the exchange z <-> 1 — z 
it is noi possible to have three linearly independent solutions with such a symmetry. When 
the series defining the hypergeometric functions do not truncate there are singularities at 
u = 4z(l — z) = 1 of the type log(l + y/1 — u). And since the reciprocal image of the 
interior of the u-unit disc consists of two disjoint lobes, in the z plane, centered around 
and 1, the change z — > 1 — z changes the determination of the hypergeometric function 
around z = 1 from what it is around z — 0. This feature can be seen explicitly for the 
cases x == 1 and x = 1/2 below. 

Before analyzing some specific solutions let us note that after the change of function 
F~(z) = (z(l — z))~*G(z) one obtains 

z 2 (z - l) 2 d 3 G + z(l - z)(K lZ + K 2 (z - l))d 2 G + (L lZ 2 + L 2 (z - l) 2 

+L 3 z(z - l))dG + (M lZ + M 2 {z - l))G = , (71) 



See 1 20 or any book on differential equations. 
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where 

K 1= K 2 = 3--, L 1 = L 2 = 1--, L 3 = 4(l-i) 2 , M 1 = M 2 = -(- - 1) . (72) 

JU JU %AJ ikj ikj 

Equations of this type appeared for correlations functions containing a level three null 
vector in minimal models (see ref. pip . There are integral representations for the solutions 

G{z) = [ dtj dt 2 t\~\t x - - z)k\'\t 2 - l)*-\t 2 - *)-(<! - t 2 )~' . (73) 

JCi Jc 2 

These integrals can however diverge, depending on the contours C\ and C 2 and the value 
of x. 



7.2 The level 1 correlators 

For x = 1, i.e. k = 1, we are able to write the solutions of in a closed form: 

F - w .^ +ffl ( 1 + _^) +92 ( 1 + M^). (74) 

We then want to construct the physical correlation functions which must be single- valued 
on the whole Riemann sphere. It is therefore enough to ensure this property at the two 
singular points and 1. It is easy to see that one can only have 

F AB (z,z) = a F A \z)FP(z) + a 1 (F A \z)F B 1 \z)+F A 1 \z)FP(z)) 

+ a 2 {Ff\z)F B 2 \z) + F A 2 \z)F B °\z)) , (75) 

where the constants are independent of the indices A, B, and 

= lift) , F?\z) = , FP(z) = \ 

^(^^(-l + glogz) , F 2 {l) (z) = ^-^(1- z + (2z-l) log z) , 
F 3 (1) (z) = ^-l + 21og^) , F?\z) = ± (-A_ + f^log(l -zj) , 
= jp^z) ( z + C 1 - ^ Ml - *)) . ^) = ^ " 2 log(l - z)) .(76) 

This form of Fab differs from the usual diagonal one one finds for minimal models. 
Here there are three so far unconstrained constants (one constant is an overall normaliza- 
tion), even after requiring monodromy invariance around the three singular points. 

We still have to require crossing symmetry among the correlators: 

F ijkl (z, z) = F ikjl {l - z ,l-z), F ijkl {z, z) = z- 1/x z~ 1/x F ilkj {l/z, l/z) . (77) 

These conditions imply that ot\ = a 2 = 0. These correlators are thus free from logarithms. 

At level one the free fields (7, ip + , 77) provide a representation [0, |] of the affine 
algebra osp(2/2)i. The corresponding correlation functions are then easily calculated by 
making use of Wick's theorem. They agree with the results obtained from the above 
calculation, and are free of logarithms. 
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We then looked at x = 1/2 (k = 3/2) correlators for which we explicitly found the 
functions . Again the physical correlators were free from logarithms. It is possible to 
show that such is the case for all non- vanishing integer values of sc -1 . The solutions F A ^ 
have the same form as those for x = 1 and x = 1/2, namely a combination of rational 
functions of z and of logz and log(l — z). We were not able to determine the connection 
matrices between functions with arguments z, 1 — z and 1/z. The limiting values for 
those obtained in [^1J depend on the way the limit is taken. But it is enough to know the 
form of these matrices to implement the crossing symmetries. In doing so we also used 
the matrices connecting the tensors I a with transposed indices, and which can be found 
in the appendix. We found that ai = «2 = 0. Note that the values of k considered here 
are exactly those of equation (flCf ). We comment on this result in the conclusion. 

From the four-point correlation functions we can obtain information about the fusion 
rules for the [0,1/2] representation. As the contribution from the conformal blocks 1 
and 2 vanishes we infer that the corresponding operators do not appear in the OPE 
0[o,i/2](^)0[o,i/2](w). The only primary field appearing in their OPE corresponds to the 
operator t , which is now the identity. 

7.3 Logarithms for generic x 

We now consider the generic case where 2x~ x is not an integer. Contrary to what happens 
in the preceding section, it does not seem possible to isolate a solution which is regular 
at both z = and z = 1. A solution regular at will have a logarithmic singularity at 
z = 1 and vice-versa. This means that there is no way to get rid of logarithms if the 
correlators are not vanishing. We expect that there should be logarithmic operators, and 
the knowledge of the four-point functions allows one to find the OPE's of the fields we 
consider, including the logarithmic ones. 

7.4 Squared logarithms 

We now turn to the third and last possible scenario where 2x~ x is an integer but x~ l is not 
an integer: x~ l = n + 1/2, n6Z. The theory of fuchsian differential equations indicates 
that there might be a solution with a log 2 z dependence. A careful analysis shows that 
this is indeed the case. And again it seems impossible to isolate a logarithm-free solution. 
However one has to modify the form of Fab- Write the solutions as 

F A °\z) = F A (z) , F%\z) = F A (z) \ogz + H A (z) , 
F ( A 2 \z) = F A (z) log 2 z + 2H A (z) \ogz + G A (z) . (78) 

Then the monodromy invariant combinations around z = are 

F AB (z,z) = a Ff(z)FP(z)+a l {Ff\z)F^\z)+F^\z)F^\z)) 

+ a 2 {Ff\z)F ( i ) (z) + F A 2 \z)F^ (z) + 2F a 1 \z)f£\z)) . (79) 

We expect such a form to be also invariant around 1 without any additional constraint 
on the aj. Again, because there is no solution free of logarithms one expects a log 2 
dependence if the physical correlators do not vanish. 
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8 Conclusion 



The current superalgebra approach to gaussian disordered systems provides a unifying 
framework for determining exponents of the models and correlation functions. The loga- 
rithmic dependence of correlation functions at critical points is traced back to indecom- 
posable representations which are common in super algebras. Chiral correlation functions 
possess a logarithmic dependence. However for values of the level k for which the theories 
have good locality properties, the physical operators do not have a logarithmic depen- 
dence. This conclusion was not apparent from the start; it came out of imposing crossing 
symmetry. It is not clear at this stage what is the significance of the logarithmic operators 
in the context of the disordered systems. If they have vanishing dimension, as is the case 
here, they might be related to the density of states. 

The correlation functions with log and log 2 dependence provide examples of correlators 



whose existence was predicted in ||16|| . The negative conformal weights are natural in this 
treatment of disorder because one is considering non- unitary theories at c = 0. It should 
be also possible to obtain similar results by considering WZW models of cosets of Kac- 
Moody superalgebras. 
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9 Appendix 

The non-vanishing components of the three invariant tensors are equal to the components 
of the three vectors: 

h = (1144) + (1234)4e7 + (1324)4e7 - (1414) + (2143)4e7 

+ (2233)16e 2 7 2 + (2323)16eV - (2413)4e7 + (3142)4e7 + (3232)16eV 

+ (3322)16e 2 7 2 -(3412)4e7- (4141) -(4231)4e 7 - (4321)4e7 + (4411) (80) 

h = (1234)4e7 - (1243)4e7 + (1324)4e7 - (1342)4e7 

- (1414) + (1441) - (2134)4e7 + (2143)4e7 + (2233)32eV + (2323)16eV 
+ (2332)16e 2 7 2 - (2413)4e7 + (2431)4e7 - (3124)4e7 + (3142)4e7 

+ (3223)16e 2 7 2 + (3232)16eV + (3322)32eV - (3412)4e7 + (3421)4e7 

+ (4114) - (4141) + (4213)4e7 - (4231)4e7 + (4312)4e7 - (4321)4e7 (81) 

h = (1234) - (1243) + (1324) - (1342) + (1423) + (1432) - (2134) + (2143) 

+ (2233)8e7 + (2314) + (2323)8e7 + (2332)8e7 - (2341) - (2413) + (2431) 

- (3124) + (3142) + (3214) + (3223)8e7 + (3232)8e7 - (3241) + (3322)8e7 

- (3412) + (3421) - (4123) - (4132) + (4213) - (4231) + (4312) - (4321) (82) 
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The matrices P and Q are given by: 



/ 1 



P 











-3 -J- 

2e7 



Q 




\ 4e7 8e7 1 

The following tensors are needed in order to implement the crossing-symmetries: 



L .4 



' ]^ £ £*2 ( £ a3+ £ a 4 )+£a 3 £c« 4 jaiCt4a3"2 (84) 



The £ ai refers to the parity of the index ctj. One has the following decomposition: 

Ji = -h, J2 = -h + h - 4e 7 / 3 , J 3 = -/ 3 (85) 
K^h-Ae^h, K 2 = h-4eyl 3 , K 3 = -I 3 . (86) 

We rewrite this compactly as J — C\ I and K = C% I where are two matrices which 
are used to implement the crossing symmetry constraints. 
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